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Abstract 

We establish existence with sharp rates of decay and distance from the Chapman- 
Enskog approximation of small-amphtude shock profiles of a class of semilinear relax- 
ation systems including discrete velocity models obtained from Boltzmann and other 
kinetic equations. Our method of analysis is based on the macro-micro decomposi- 
tion introduced by Liu and Yu for the study of Boltzmann profiles, but applied to the 
stationary rather than the time-evolutionary equations. This yields a simple proof by 
contraction mapping in weighted H'^ spaces. 
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9 Apphcation to spectral stability 

1 Introduction 

We consider the problem of existence of relaxation profiles 

(1.1) U{x,t) = tj{x-st), lim tJ{z) = U± 

2— >±00 

of a semilinear relaxation system 

(1.2) Ut + F{U)^ = Q{U), 
in one spatial dimension, with the following structure: 

Assumption 1.1. (HI) The flux F is linear in U , so that F{U) = AU for some constant 
matrix A. 

(H2) There are linear coordinates such that 

(1.3) ^^(:). Q^(» 

(H3) q has nondegenerate equilibria parametrized byu; more precisely there are a smooth 
function from R" to R'' and 9 > such that for all u: 

(1.4) q(u,v^,{u)) = 0, Re a{dyq{u,V:t{u))) < —9, 
(t(-) denoting spectrum and Re the real part. 

Common examples are discrete kinetic models obtained by discrete velocity or other 
approximation from continuous kinetic models such as Boltzmann or Vlasov-Poisson equa- 
tions; for example, Broadwell and other lattice gas models [PI]. Other examples are the 
semilinear relaxation schemes introduced by Jin-Xin [JX] and Natalini [N] for the purpose 
of numerical approximation of hyperbolic systems. Here, we are thinking particularly of the 
case n bounded and r >> 1 arising through discretization of the Boltzmann equation, or 
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the case r — > oo arising in Boltzmann itself; that is, we seek estimates and proof independent 
of the dimension of v. 

For fixed n, r, the existence problem has been treated in [YZllMaZT] under the additional 
assumption 

(1.5) det{dF - si) ^ 

corresponding to nondegeneracy of the traveling-wave ODE. However, as pointed out in 
|MaZ2t IMaZ3] . this assumption is unrealistic for large models, and in particular is not 
satisfied for the Boltzmann equations, for which the eigenvalues of dF are constant particle 
speeds of all values. Our goal here, therefore, is to revisit the existence problem without the 
assumption (jl.Sp . with the eventual aim being to establish a simple proof of existence of 
small-amplitude Boltzmann profiles. Of course, existence of such was established some time 
ago in |CN] : however, the proof is rather complicated, involving detailed resolvent estimates 
in weighted L°° spaces in spatial and velocity variables, and so it seems of use to seek a 
simpler approach based on weighted spaces and standard energy estimates. 

Our method of analysis is motivated by the "macro-micro decomposition" technique 
introduced by Liu and Yu [LYj . in which fluid (macroscopic, or equilibrium) and transient 
(microscopic) effects are separated and estimated by different techniques. This was used 
in |LYj to show by a study of the time-evolutionary equations that the Boltzmann pro- 
file constructed in |CNj has nonnegative probability density, that is, to show positivity of 
Boltzmann profiles assuming that such a profile exists. 

Our approach here is very much in the spirit of that of [LYJ, based on approximate 
Chapman-Enskog expansion combined with Kawashima type estimates (the macro-micro 
decomposition of the reference), but carried out for the stationary (traveling- wave) rather 
than the time-evolutionary equations, and estimating the finite-dimensional fluid part using 
sharp ODE estimates in place of the energy estimates of |LYj. In this latter part, we are 
much aided by the more favorable properties of the stationary fluid equations, a rather 
standard boundary value ODE system, as compared to the time-evolutionary equations, a 
hyperbolic-parabolic system of PDE. 

Our main result is to show existence with sharp rates of decay and distance from the 
Chapman-Enskog approximation of small- amplitude quasilinear relaxation shocks in the 
general case that the profile ODE may become degenerate. See Sections [2] and [3] for model 
assumptions and description of the Chapman-Enskog approximation, and Section [5] for 
a statement of the main theorem. In the present, semilinear case, a simple contraction- 
mapping argument suffices; the quasilinear case is treated by Nash-Moser iteration in 
|MeZl| . In |MeZ2] . we show that the argument of this paper carries over with minor 
modifications to the infinite-dimensional Boltzmann equation with hard potential to yield 
existence of small-amplitude Boltzmann shock profiles, recovering and slightly sharpening 
the results of [CNj . This in a sense completes the analysis of plY] . providing by a common 
set of techniques both existence (through the present argument) and (through the argu- 
ment of [LYJ ) positivity. At the same time it gives a truly elementary proof of existence of 
Boltzmann profiles. 
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Finally, we note that spectral stability has been shown for general small-amplitude quasi- 
linear relaxation profiles in |MaZ3j . without the assumption (jl.5p . under the assumption 
that the profile exist and satisfy exponential bounds like those of the viscous case. The 
results obtained here verify that assumption, completing the analysis of [MaZ3] . It would 
be very interesting to continue along the same lines to obtain a complete nonlinear stability 
result as in [MaZlj . in particular for Boltzmann shocks. 

Existence results in the absence of condition (II. Sp have been obtained in special cases 
in [MaZ5l IDYj by quite different methods (for example, center-manifold expansion near 
an assumed single degenerate point [DY]). However, the decay bounds as stated, though 
exponential, are not sufficiently sharp with respect to e for the needs of |MaZ3j . More 
important, the techniques used in these analyses do not appear to generalize to the infinite- 
dimensional (e.g., Boltzmann) case. 

2 Model, assumptions, and the reduced system 

Taking without loss of generality s = 0, we study the traveling- wave ODE 
(2.1) AU' = Q{U), 



(2.2) [/=('"), A = constant, Q = ( . 

governing solutions of (jl.lh . where q satisfies (II. 4p . 
We use the notations 



(2.3) A 



'All Ai2 

^A2l A22^ 

(2.4) f{u,v) := Aiiu + Ai2V. 



We make the standard assumption of simultaneous symmetrizability [Y]: 

Assumption 2.1. (SS) There exists a smooth, symmetric and uniformly positive definite 
matrix S{U) such that 

i) for all U , S(U)A is symmetric, 

a) for all equilibria = {u,v^:{u)), SdQiU^,) is symmetric nonpositive with 

(2.5) dimker = dimkerdQ = n. 

We also make the Kawashima assumption of genuine coupling [K]: 

Assumption 2.2. (GC) For all equilibria = (u, f*(u)), there exists no eigenvector of 
A in the kernel of dQ(U^). Equivalently, given Assumption \2.1\ (see JK^), there exists in a 
neighborhood J\f of the equilibrium manifold a skew symmetric K = K{U) such that 

(2.6) Re {KA - SdQ){U) >9>0 
for all U eM. 
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Recall from [Y] (see aslo Section [3|), that the reduced, Navier-Stokes type equations 
obtained by Chapman-Enskog expansions are 

(2.7) Mu)' = {K{u)uy, 
where 

(2.8) /*(n) := f{u,v^{u)) = Auu + Ai2V^{u), 

(2.9) h{u) := -Ai2c4u) 

with 

c*(u) := di,q~^{u,v^{u)) 

(2.10) / N 

{A21 + A22dv^{u) - dv^:{u){Aii + Ai2dv^{u)]. 

Note also, by the Implicit Function Theorem, that 

dv^,{u) = —dyq^^duq{u,v^,{u)). 

For the reduced system (j2.7p . simultaneous symmetrizability becomes: 

(ss) There exists s{u) symmetric positive definite such that sdf^, is symmetric and sb^, 
is symmetric positive semidefinite. 

We have likewise a notion of genuine coupling [K] : 

(gc) There is no eigenvector of d/* in ker 6^, . 

We note first the following important observation of [Y] . 

Proposition 2.3 (H)- -^^^ (|2.1|) as described above be a symmetrizable system satisfying 
the genuine coupling condition (GC). Then, the reduced system (j2.7p is a symmetrizable 
system satisfying genuine coupling condition (gc). 

Proof. We only give here a sketch, mentioning the key points and refereeing e.g. to [MaZ3] 
for details. Fixing U, one is reduced to constant matrices and linear algebra, with matrices 



^A2i A22J ' VQ21 Q22 

and symmetrizer S. With 

P=(^fr V = -Q2iQ2 



V Id) ' ~^22 V21, 

the change of unknowns U = PU transforms the problem to an equivalent one with matrices 
A = P-^AP, Q = P"^QP and symmetrizer S = P*SP, with 



Q 





Q22 
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Therefore S is block diagonal (52i = and 5i2 = 0), SuAu is symmetric, SuA^ = 
(522^21 )* and S22Q22 is definite negative. Next, the associated matrix b is : 

6= -Ii2Q22^l21- 

Thus 5*11 is definite positive, symmetrizes An and 

Sub = —5'iij4i2 (522^22^ ) ^'S'22^21 = —(522^21)* ('S'22Q22^)~^'^22^21 

is symmetric and nonnegative. Noticing that 

(2.11) An =^11 + ^12^, 

(2.12) b = -A^2Q22 (^21 + ^22^ - y{Au + ^12^)) = 6* 

this implies that the property (ss) is satisfied with symmetrizer s = Su (In terms of the 
original matrices, s = Sn + S12V = Su + S12V + V*{S2i + 5'22^)) since ^21 + ^22^ = as 
a consequence of the block diagonal structure of S). 

Similarly, the property (GC) is transported to the system (A,Q), meaning that 

Auu = Xu, A21U = =^ u = 0. 

The symmetry property of Sub implies that 

ker b = ker A21 

and the property (gc) immediately follows. □ 

Besides the basic properties guaranteed by Proposition 12. 3| we assume that the reduced 
system satisfy the following important additional conditions. 

Assumption 2.4. (i) The matrix 6*(it) has constant left kernel. 

(a) For all values of u, ]<ienr^:df^{u) Piker 6* (tt) = {0}, where 7r*(u) is the zero eigenpro- 
jection associated with 6,,(n). 

The importance of Assumption 12.41 in the present situation is that it ensures that the 
zero-speed profile problem for the reduced system, 

(2.13) f*{u)' = {b*{u)u')' , lim u{z) = u± 

2:— >±oo 

or, after integration from —00 to x, 

(2.14) 6,(u)n' = /,(n) -/,(n±), 

may be expressed as a nondegenerate ODE in U2, coordinatizing u = (ui, n2) with ui = vr^u 
and U2 = {I — 7r*)n. 

Condition (i) was pointed out in [Zej, condition (ii) in [MaZSl [ZT j IGMWZ] . 
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Remark 2.5. Assumption 12.41 is also central to the linearized stability analysis of general 
Navier-Stokes type equations in \Z2\ IGMWZj . It appears to be independent from the 
genuine coupling conditions (GC), (gc), except in the special case that u is scalar, for which 
(GC), (gc) reduce to kerb,, = 0. It is satisfied for the important example of Boltzmann 
equations, which is our main motivation. 

Next, we assume that the classical theory of weak shocks can be applied to (|2.13p . 
assuming that the flux /* has a genuinely nonlinear eigenvalue near 0: 

Assumption 2.6. In a neighborhood U^, of a given base state uq, df^, has a simple eigenvalue 
a near zero, with a{uo) = 0, and such that the associated hyperbolic characteristic field is 
genuinely nonlinear, i.e., after a choice of orientation, Va • r{uo) < 0, where r denotes the 
eigendirection associated with a. 

Remark 2.7. Assumption [2]6] is standard, and is satisfied in particular for the compressible 
Navier-Stokes equations resulting from Chapman-Enskog approximation of the Boltzmann 
equation. Assumptions 12.11 and [2] are verified in for a wide variety of discrete kinetic 
models. Assumptions [274] and \2M on the reduced equations must be checked in individual 
cases. 

3 Chapman— Enskog approximation 

Integrating the first equation of (12. Ih and noticing that the end states {u±,v±) must be 
equilibria and thus satisfy v± = v^{u±), we obtain 



The idea of Chapman-Enskog approximation is that v — v^{u) is small (compared to the 
fluctuations u — u±). Taylor expanding the second equation, we obtain 



(3.1) 



Aiiu + Ai2v = f*{u±), 
A21U' + A22V' = q{u, v). 



Because / is linear, the first equation reads 



(3.2) 



/*(n) + Ai2{v - v^{u)) = /*(u±). 



{A21 + A22dv^.{u))u + A22{v - v^{u)y = dvq{u,v^{u)){v - w*(n)) 

+ 0{\v-v4u)\''), 



or inverting d^q 



(3.3) 




The derivative of (|3.2p implies that 



[Ann + Ai2dv^{u))u' = 0{\{v - f*(n))'|). 
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Therefore, ()3.3p can be replaced by 

(3.4) V - v^u) = c,{u)u' + 0{\v - 7;,(u)|2) + 0{\{v - v,{u))'\), 



where is defined at (j2.10p . Substituting in (j3.2p . we thus obtain the approximate viscous 
profile ODE 

(3.5) h{u)u' = Mu) - Mu±) + 0{\v - v,{u)\^) + 0{\{v - t;*(u))'|), 
where 5* is as defined in (|2.9p . 

Remark 3.1. The above calculation leaves a great deal of flexibility in the choice of 6^, 
satisfying (|3.5p . namely it is only specified modulo multiples of 

All + Ai2dv^,{u), 

as we used when passing from (|3.3p to (|3.4p . However, we have chosen to use the standard 
definition ()2.8p of 6* because it is the natural choice which is invariant by change of variables 
(see (j2.12p ) and it is known by Proposition 12.31 and by the explicit example of Boltzmann to 
have good properties. But, it might be, for example, that a different representative could 
be strictly parabolic, so slightly easier to handle. This seems to be just a curiosity, as 
the analysis is already sufficient to treat the standard case. But, it is interesting from the 
viewpoint of the Chapman-Enskog expansion and possible alternative representations. 

Motivated by (j3.3p - (j3.5p . we define an approximate solution (uns^^ns) of (|3.ip by 
choosing n^r^ as a solution of 

(3.6) b:,{uNs)uNs = f*iuNs) - f*iu±), 
and vns as the first approximation given by (j3.3p 

(3.7) VMS - v^{uNs) = c*{uns)u'ns- 

Small amplitude shock profiles solutions of (j3.6p are constructed using the center man- 
ifold analysis of [Pej under conditions (i)-(ii) of Assumption 12. 4t see discussion in |MaZ5j . 

Proposition 3.2. Under Assumptions \2. ^ and \2.4\ in a neighborhood of {uo,uo) in x 
R", there is a smooth manifold S of dimension n passing through {uo,uq), such that for 
(n_,ii+) E S with amplitude e := — n_| > sufficiently small, and direction {u+ — u-)/e 
sufficiently close to r(no), the zero speed shock profile equation (j3.6p has a unique (up to 
translation) solution u^s The shock profile is necessarily of Lax type.' i.e., with 

dimensions of the unstable subspace of df^{u-) and the stable subspace of df^,{uj^) summing 
to one plus the dimension of u, that is n + 1. 

Moreover, there is 6 > and for all k there is independent of (u_,u+) and e, such 
that 

(3.8) \d',{uNS-u±)\<Cke''+'e-''\''\, x^O. 
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We denote by 5+ the set of m+) G S with amphtude e := — U-\ > sufficiently 
small and direction {u+ — n_)/e sufficiently close to r(no) such that the proffie uns exists. 
Given (n_,M+) € 5+ with associated profile uns^ we define v^s by (|3.7|) and 

(3.9) Uns ■= iuNS,VNs)- 

It is an approximate solution of (jS.ip in the following sense: 
Corollary 3.3. For (n_,M+) G 5+, 

(3.10) ^iiniV5 + AuVNS - f*{u±) = 
and 

satisfies 

(3.11) \d^T^v(.x)\ < Cfce'^+^e-^'l^l, x^O 
where Ck is independent of {u-,u^) and e = — n_|. 

Proof. Given the choice of vns^ the first equation is a rewriting of the profile equation (|3.6p . 
Next, note that 

VNS - v*(uNs) = 0{\u'ns\), {vns - v^iuNs))' = 0{\u%s\) + Od-u^sP), 

where here O(-) denote smooth functions of uns and its derivatives, which vanish as indi- 
cated. With similar notations, the Taylor expansion of q and the definition of vns show 
that 

7^^, =0{\VNS - V*{UNS)\'^) + 0{\{VNS " t^* (""Afs))'! ) 

+ dv^{uNs {All + Ai2dv^{uNs))u'Ns- 

Moreover, 

(^11 + Ai2dv^{uNs))uNs = {f*iuNs)y = {b*{uNs)uNs)y 

= 0{\u'^s\')+0{\u%s\)- 

This implies that 

n, = 0{\u'^s\') + 0{\u%s\)- 
satisfies the estimates stated in ()3.1ip . □ 

Remark 3.4. One may check that if we did not include the correction from equilibrium 
on the righthand side of (|3.7p . taking instead the simpler prescription v^s = v^{uns) as 
in [LY] . then the residual error that would result in ()3.10p would be too large for our later 
iteration scheme to close. This is a crucial difi^erence between our analysis and the analysis 

of (Hj. 
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4 Statement of the main theorem 

We are now ready to state the main result. Define a base state Uq = {uo,v^{uq)) and a 
neighborhood U = x V. 

Theorem 4.1. Let Assumptions (SS), (GC), and \2.4\ hold on the neighborhood U of Uq, 
with Q £ C°° . Then, there are Eq > and 6 > such that for {u-,u^) £ S+ with amplitude 
£ := u_| < eo, the standing-wave equation (12.11) has a solution U inU, with associated 

Lax-type equilibrium shock (n_,n+), satisfying for all k: 

\d^AU -UNs)\<Cke''^^e-''\^\, 
(4.1) \dt{u - n±)| < Cfce^+^e-'^^l^l, x ^ 0, 

\dl{v-v,{u)\ < Cfce'^+V'^^l^l, 

where tJNS = {uns,'(^ns) is the approximating Chapman-Enskog profile defined in (|3.9p . 
and Ck is independent of e. Moreover, up to translation, this solution is unique within a 
ball of radius ce about Uns in norm \\-\\l2 -^e~^\\dx •\\l2 +£~'^\\d^ for c > sufficiently 
small. (For comparison, Uns ~ U± is order e^^'^ in this norm, by l\4:.l\\ (ii)-(iii).) 

Bounds (|4.ip show that (i) the behavior of profiles is indeed well-described by the Navier- 
Stokes approximation, and (ii) profiles indeed satisfy the exponential decay rates required 
for the proof of spectral stability in |MaZ3j . Prom the second observation, we obtain 
immediately from the results of |MaZ3j the following stability result, partially generalizing 
that of |LY] for the Boltzmann equations 

Corollary 4.2 ( [MaZSj ). Under the assumptions of Theorem \4-ll the resulting profiles 
U are spectrally stable for amplitude e sufficiently small, in the sense that the linearized 
operator L := dxA{U) — dQiU) about U has no L^ eigenvalues A with Re A > and A 7^ 0. 

The remainder of the paper is devoted to the proof of Theorem 14.11 

Remark 4.3. Theorem 14. 1 1 yields uniqueness only among solutions close to the Chapman- 
Enskog approximant Uns- The stability result of Liu-Yu |LY] should give uniqueness 
among solutions in a ball of small but 0(1) radius, assuming that they have zero relative 
mass compared to Uns- Indeed, it should be possible to upgrade this to general-mass 
perturbations to obtain ultimately a full 0(1) uniqueness result. Stability with respect to 
general-mass perturbations is an important open problem. 

^ Liu and Yu prove the stronger result of linearized stability with respect to zero- mass perturbations that 
are sufficiently small in an appropriate norm. 
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5 Outline of the proof 



5.1 Nonlinear perturbation equations 

Defining the perturbation variable U := JJ — Uns, and expanding about Uns, we obtain 
from (j3.ip the nonlinear perturbation equations 

(5.1) Aiiu + Ai2V = 

(5.2) A21U' + A22V' - dqiUNs)U = -n^ + N{U) 

where the remainder N(U) is a smooth function of Uns and U, vanishing at second order 
at [/ = 0: 

(5.3) N{U)=Af{UNS,U) = 0{\U\''). 
We push the reduction a little further, using that 

(5.4) M := dq{uNs,VNs) - dq{uNS,v*{uNs)) = 0{\vns - v^{uns)\)- 
Therefore the equation reads 



(5.5) 





-7^^, + MU + N{U) 



where 

(5.6) Q21 = duqiuNS,v*{uNs)), Q22 = d^q{uNS,v*iuNs))- 
Differentiating the first line, it implies that 

(5.7) LIU := AU' - dQ{uNS,v.{uNs))U = (^_^^ ^ ^ ^^^^ 

The linearized operator Adx — dQ{U) about an exact solution U of the profile equations 
has kernel U', by translation invariance, so is not invertible. Thus, the linear operators 
LI and CI are not expected to be invertible, and we shall see later that they are not. 
Nonetheless, one can check that is surjective in Sobolev spaces and define a right inverse 
£1{CI)^ = I, or solution operator (Cl)^ of the equation 

(5.8) ClU = F-={^g 

as recorded by Proposition 15.21 below. Note that is not surjective because the first 
equation requires a zero mass condition on the source term. This is why we solve the 
integrated equation ()5.5p and not (|5.7p . 
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To define the partial inverse we specify one solution of (jS.Sp by adding the co- 

dimension one internal condition: 

(5.9) 4 • m(0) = 
where 4 is a certain unit vector to be specified below. 

Remark 5.1. There is a large flexibility in the choice of 4. Conditions like (j5.9p are 
known to fix the indeterminacy in the resolution of the linearized profile equation from 
(j3.6|) and it remains well adapted in the present context, see section [7] below. A possible 
choice, would be to choose 4 independent of e and parallel to the left eigenvector of df^{uo) 
for the eigenvalue (see Assumption I2.6p . 

5.2 Fixed-point iteration scheme 

The coefficients and the error term TZv are smooth functions of uns and its derivative, thus 
behave like smooth functions of ex. Thus, it is natural to solve the equations in spaces 
which reflect this scaling. We do not introduce explicitly the change of variables x = ex, 
but introduce norms which correspond to the usual norms in the x variable : 

(5.10) 11/11^. =e5||/||^2+e-^||a,/||i2 + ... + e5-||9^/||^2. 

We also introduce weighted spaces and norms, which encounter for the exponential decay 
of the source and solution: introduce the notations. 

(5.11) < X >:= (x^ + 1)1/2 

For 5 > (sufficiently small), we denote by ^ the space of functions / such that e^^^^^ f G 
H'^ equipped with the norm 

(5.12) =^^E^~'II^''^'^^'/II^- 

k<s 

Note that for 6 < 1, this norm is equivalent, with constants independent of e and 6, to the 
norm 

Proposition 5.2. Under the assumptions of Theorem \4-l\ there are constants C, eo > 
and 5o > and for all e £]0, e^], there is a unit vector ii, such that for e g]0, Eq], S G [0, 6o], 
f G s> 9 ^ 5 operator equations (15. 8|) (15. 9p has a unique solution U G ^, 
denoted by U = {Cl)^F, which satisfies 

(5.13) <Ce-W\f\\HS +\\g\\^, ). 

Moreover, for s > 3, there is a constant Cg such that for e g]0, Eq] ond / G H^^^ , 
g G H^^ the solution U = {C^^F G H^^ and 

(5.14) \\{Cl)^F\\^ < Ce-W\f\\^sy + ll^ll^. ) + Cs\\{CtyF\\ 
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The proof of this proposition comprises most of the work of the paper. Once it is 
estabhshed, existence follows by a straightforward application of the Contraction-Mapping 
Theorem. Defining 

^^•^^^ (-7^, + ML/ + iV(C/))) ' 

we reduce (j5.7p to the fixed-point equation 

(5.16) TU := U. 

5.3 Proof of the main theorem 

Proof of Theorem \4-l\ The profile ums exists if e is small enough. The estimates (|3.8p 
imply that 

(5.17) \\UNS - '"±11^1,5 < CsE 

with Cs independent of e and 5, provided that 5 < 6/2. Similarly, (jS.lip implies that 

(5.18) \\n4Hi,<Cse\ 
and (|5.4p implies that 

(5.19) \\M\\hi^ <Cse'^. 

Moreover, with the choice of norms ()5.10p . the Sobolev inequality reads 

(5.20) II-uIIlcx, < C\\u\\hi < C||n||^i 

with C independent of e. Moreover, for smooth functions <I>, there are nonlinear estimates 

(5.21) \mu)\\Hi < C{\\u\\loo) Mhi- 
which also extend to weighted spaces, for 5 < 1: 

(5.22) mu)\\Hi^^<C{\\u\\L^) \\u\\hi,. 

In particular, this implies that for s > 1, 5 < min{ 1,0/2} and e small enough: 
I|MC/||h|,,<C(||M||^iJ|C/||h;,, + ||M||^,.J|[/||^.^ 
<e'{C\\U\\Hi^ + Cs\\U\\Hi^) 
where the first constant C is independent of s. Similarly, 

(5.24) mmiHi, < C{\\U\\l^)\\U\\hiJU\\hi,. 
Combining these estimates, we find that 

\\TU\\h^^ < e-\Cse^ + Ce^UWhs + + C||C/||^i ||C/||j^. ), 

^ £,0 £,0 £,0' 

that is 

(5.25) W'^UWhIs ^ + + ^^'IMhI^WWhIs + Css\\U\\h^^^ 
provided that £ < £o, ^ ^ min{l,6'/2} and ||C/||ioo < 1. 
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Consider first the case s = 2. Then, T maps the ball 8^5 = {\\U\\jj2 < to itself, 

if e < El where £i > is small enough. Similarly, 

(5.26) \\TU-TV\\Hi^<Ce-\e^ + \\U\\Hi + \\V\\Hi)\\U-V\\Hi^, 

provided that ||?7||/,oo < 1 and ||l/||j;,oo < 1, from which we readily find that, for e > 
sufficiently small, T is contractive on B^^s, whence, by the Contraction-Mapping Theorem, 
there exists a unique solution of (j5.16p in B^^g for e sufficiently small. 
Moreover, from the contraction property 

\\u' - T(o)|U| = wrm - TmiHi < cWu'Whi, 

with c < 1, we obtain as usual that IIC/^'^^H 71-2 < C||T(0)||Lr2 , whence 

£,S £,S 

(5.27) \\U''\\h2^<C£'^. 

by (|5.25p . In particular, = O(e^) in and by the Sobolev embedding 

(5.28) ||e^^<^)c/^||ioo =0(e2), \\e^^^^'^ d^U^U-o = 0{e^). 

For s > 3, the estimates (15.25P show that for e < ei independent of s, the iterates T"'(0) 
are bounded in ^, and similarly that T"'(0) — T(0) = 0{e'^) in implying that the 
limit U belongs to ^ with norm 0{e'^). Together with the Sobolev inequality (j5.20p . this 
implies the pointwise estimates (14. ip . 

Finally, the assertion about uniqueness follows by uniqueness in Bc^^s for the choice 5 = 
and c > sufficiently small (noting by our argument that also Bce^s is mapped to itself for 
e sufficienty small, for any c > 0), together with the observation that phase condition ()5.9p 
may be achieved for any solution U = Uns + U with 

||f/||L- < ce^ « U'^si^) ~ 

by translation in x, yielding Ua{x) := U{x + a) = Uns{x) + Ua{x) with 

Ua{x) := Uns{x + a) - Uns{x) + U{x + a) 

so that Ua{0) ~ (a + o{l))IJ'j^g{0) and £s • Ua{0) ~ ■ u'j^g{0), which may be set to zero by 
appropriate choice of a, by the property is ■ u'j^si^) ^ following from our choice of ie (see 
Remark [511. □ 



It remains to prove existence of the linearized solution operator and the linearized bounds 
()5.14p . which tasks will be the work of the rest of the paper. We concentrate first on 
estimates, and prove the existence next, using a viscosity method. 
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6 Internal and high frequency estimates 



We begin by establishing a priori estimates on solutions of the equation (jS.Sp This will 
be done in two stages. In the first stage, carried out in this section, we establish energy 
estimates showing that "microscopic" , or "internal" , variables consisting of v and derivatives 
of {u, v) are controlled by and small with respect to the "macroscopic" , or "fluid" variable, 
u. In the second stage, carried out in Section [71 we estimate the macroscopic variable u by 
Chapman-Enskog approximation combined with finite-dimensional ODE techniques such 
as have been used in the study of fluid-dynamical shocks; see, for example, |MaZ41 IMaZ5[ 
[ZHIZ2 1IUMWZ] . 

6.1 The basic estimate 

We consider the equation 

A11U + A12V \ (f 



f6 'U 

* ' \A21u' + A22V' - dq{uNS,v*{uNs))U J \g 

and its differentiated form: 

7' 



(6.2) AU' - dQiuNs, V, {uNs))U = 



The internal variables are U' = {u',v') and v where 

(6.3) v:=v+pu, p = d^q~^duq{uNS,'V:t:{uNs)) = -dv^{uNs) 

is the linearized version of -u — (n) . 



Proposition 6.1. Under the assumptions of Theorem \4-l\ for there are constants C, Eq > 
and 60 > such that forO < e < Eq andO < S < 60, f & ^eS' 9 ^ ^1 5 '^'^dU = {u,v) & s 
of (j6.ip satisfies 



(6.4) \\U'\\^, +\\v\\^, <C{\\{f,f'J",g,g')h2+e\\u\\^,). 

e.6 £,S -•'J £,S 

Making the change of variables v ^ v, or U ^ U = {u,v) and denoting U = P{uns)U, 
we obtain an ODE 



(6.5) AU' -QU = F + CU, 
where 

(6.6) A = P-'AP, Q = p-'dQP = (° ~° J 
with Q22 = dyq{uNS,'V:t:{uNs), 



9 + R2if' 
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with R21 = d^q ^duq{uNS,v*{uNs), and 



(6.8) C = -P-'AP' = ( ll j = Oiu'^s) = e'C. 



f' + eh 
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The equation (|6.5p reads 

(6.9) AU' -QU = F 
with 

(6.10) h = eCnu, ~g = 9 + i?2i/' + e^Csm. 

We first prove the estimate ()6.4p for 5 = 0. Dropping hats and tildes, the ODE reads 

(6.11) AU'-QU = F^ Q=f^ y + 



vO Q227' V a 

The matrix ^4 = A{x) has end points values A± at ±00 and satisfies estimates 
(6.12) \d^M-A^)\<Cke'' 

with Cfc independent of e. There are similar estimates for Q22- Moreover, A and Q are 
simultaneously symmetrizable by some S = S{uns), since this property is unaffected by 
coordinate changes. S is necessarily block-diagonal, SA and 



SQ 




q 



are symmetric with q negative definite. Likewise, the genuine coupling condition still holds, 
which, by the results of [K], is equivalent to the Kawashima condition, and there is a smooth 
K = K{ums) = —K* such that Re {KA — SQ) is definite positive. Therefore, there is c > 
such that for all e < eo and x G R: 

(6.13) q < -eld, Re [KA - SQ) > cid. 

Lemma 6.2. There is a constant C such that for e sufficiently small, f G H^, g G H^, 
he andU £ satisfying ([6lT]) . one has 

(6.14) \\U'h2 + \\vU2 < C{\\f\\H2 + \\h\\H. + \\g\\H. + £11^x1^2). 
Proof. Introduce the symmetrizer 

(6.15) S = dloS + d^oK-XS. 
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One has 



Re 52 o 5 o (Ad^ -Q) = ^d^ o (SA)' o - o SQ o - Re 8^ o (SQ)' 
Re o K{Ad^ - Q) = oRe KAo - Red^ o KQ 
Re SiAd^-Q) = ^{SAy - SQ. 



Thus 



Re 5 o {Ad^ - Q) =d^ o (Re AK - SQ) o + XSQ 



+ ^d^ o (SA)' o - ^X{SAy -Red^o {SQ)' - Re o KQ. 



Therefore, for U G H'^{R), (|6T3]) imphes that 

Re {SF,U)l2 > c\\d^U \\l2 + Xc\\v\\l2 
1 



Taking 



^\\{SA)'U^{\\dMl, + X\\U\\l,) 

\\{SQy\\L^\\U\\L2\\d^U\\L2 - \\K\\L^\\d^U\\L2\\qv\\L2. 



2 2 

~ L°° \\q\\l° 

c 



and using that 

(6.16) + II(SQ)'||l- = 0{e^) 

yields 

llC^'lli^ + ll^^lli^ < Re {SF, U)l2 + em\U\\l, + \\U'\\l,). 
In the opposite direction, using the block structure of S, 

Re {SF,U)l2 <\\d^U\\L2{\\d^{SF)\\L2 + \\K\\Lo.\\F\\L2) 

+ X{e\\Su\\L^\\u\\L2\\hh2 + WiSnuYWMlL^ 



+ ll'S'22||L°°|b||L2||9llL2)- 

Using again that the derivatives of the coefficients are O(e^), this implies that 

Re {SF,U)l2 <{\\f\\m + \\h\\m + \\g\\m)\\U'\\L2 

+ e|l^llL2|kllL2 +e^ll/llL2||n||i2 + lbllL2||t;||L2, 

The estimate (|6.14p follows provided that e is small enough. 

This proves the lemma under the additional assumption that U G H^. When U ^ H^, 
the estimates follows using Priedrichs mollifiers. □ 
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Proof of Proposition \6.1\ Consider the system (|6.9p Because the coefficients are functions 
of uns and its derivatives, there holds 

11^5 < C'e||n'||/^2 +e||n||j^2 

llff,5'llL2 <C(||/',r,<7,5'llL2+e'h'llL0 

and 

Therefore, the bounds (j6.14p for {U',v) imply that 

(6.17) \\U'\\^, + \\v\\^, < C{\\{fJ',f",h,h',g,g')\\L2+e\\u\\^,). 

Multiplying by £2 and using the estimates of h and g above, yields ()6.4p for 5 = 0. 
For 6 > small, consider U"^ = e^^^^'^U. Then, satisfies 



(6.18) ClU" 



w I ' 



with = e'^'^<^>/ and g"" = e^''<^>c/ + e6{xy{A2iu'" + ^22?^'")- We note that, 
\\U'h2 < m^Yhl+eWUniLp, \\vh\ < II^'IIlI, 

£,0 ^ ^ £,0 £ 

\\r,{n',in"\\^<\\u;f'j")\\L\, 

\\9'"A9n'\\Ll<\\i9,9')\\L-^+e6\\{U,U')h. 

We use the estimate (|6.4p with 5 = for , and the Proposition follows provided that 5 
is small enough. □ 

6.2 Higher order estimates 

Proposition 6.3. There are constants C , Eq > 0, 5q > and for all k > 2, there is Ck, 
such that < e < Eq, 5 < 5q, U G g, f e H^~^^ and g G g satisfying (|6.1ip satisfies: 

||a,^C/'||^2 +119,^^1^2 < C\mfJ'J",9,9')\\L-^ 

(6.19) ' 

+ e^Ck{\\U'\\rrk-i +e\\v\\jrk-i +e\\u\\L2 ) 

£,5 e,5 

Proof. Differentiating (|6.ip k times, yields 

(6.20) Ad^U'' - dQ{uNs, v4uNs))dccU^ = (^gkf^^^ > 
where 

fe-i 
1=0 
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Here we have used that dq{uNS:V^{uNs)U = Q22V- The estimate yields 
\\d^U'\\L\ + \\dlv + pdlu\\L2^ < C{\\d^,{fj'j\g,g')\\L2 

£,6 £,6 ^ £,6 

+4d^u\\^2 + \\d^rk\\L'2 + ||rfc||i2 ), 

for < k < s, with ro = when k = 0. Since Q is a function of uns^ its k — l-th derivative 
is 0(e'^^'+^) when k — I > 0. Therefore: 



Similarly, for k = 1 



IdxTkh^ ^ + lkfc|lL2, < Ck£ (\\v'\\fjk-i +e||?)||i2 ). 



IdxVkWi'^ ^ < \\dxv + pdxu\\^2 + Ce ||u||^2 

£,o £,d £,d 



and for k > 2: 



\d%\\L2 < \\d^v+pd^u\\L2 +Ck(,e''\\u'\\j^,-2+e''+^\\u\\L2 ). 

£,6 £,d £.6 



□ 



7 Linearized Chapman— Enskog estimate 
7.1 The approximate equations 

It remains only to estimate ||u||r2 in order to close the estimates and establish ()6.4p . To 
this end, we work with the first equation in (j6.ip and estimate it by comparison with the 
Chapman-Enskog approximation (see the computations Section [3]) . 
From the second equation 

^21?^' + ^221^' - g = duQU + d^qv = dvQV, 
where we use the notations v of Proposition 16.1^ we find 

(7.1) V = dyq'^(^{A2i + A22dydv^{uNs))u + M2v' - g^. 
Introducing v in the first equation, yields 

(^11 + Ai2dv^{uNs))u + Ai2V = /, 

thus 

{All + Ai2dv^{uNs))u' = f - Ai2v' - d'^v^{uNs){u'^s^'^)- 
Therefore, ()7.ip can be modified to 

(7.2) V = c^{uNs)u' + r 
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with 

r = d~'^q{uNS,v*{uNs)) (^22('v)' - 9 

This imphes that u satisfies the hnearized profile equation 

(7.3) Ku' - df,u = Ai2r - f 
where h = hiuNs) and df^ := d/*(-UArs) = ^ii + A^dv^iuNs)- 
7.2 estimates and proof of the main estimates 

Proposition 7.1. The operator b^dx — df^ has a right inverse {b^dx — df*)^ satisfying 

(7.4) \\{Kdx-df,)^h\\^2 <Ce~'\\h\\^2 

uniquely specified by the property that the solution u = {b^^^dx — df'^)^h satisfies 

(7.5) 4 • u{0) = 0. 
for certain unit vector i^. 

Taking this proposition for granted, we finish the proof of the main estimates in Propo- 
sition 15. 2[ 

Proposition 7.2. There are constants C, Eq > and 6o > such that for e G]0,eo], 
5 G [0,6o], f G H^g, g E H^^ and U £ H^^ satisfying dEHD and I^M) 

(7-6) \\U\\^, <Ce-\\\f\\H.^ + \\g\\^, ). 

Proof. Going back now to (|7.3p . u satisfies 

Ku' - df,u = 0{\v'\ + \g\ + l/'l + eM) " /, 
If in addition u satisfies the condition (j7.5p then 

(7.7) \\u\y < Ce-\\\v'\\^. + \\{fj',g)\\^. +e^u\y ). 

£ ,0 £,0 £,0 £,0 ' 

By Proposition 16.11 and Proposition 16.31 for /c = 1, we have 
(7-8) +\\v\\^, <C{\\{fj\f",g,g')\\L2+e\\u\\^,). 
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Combining these estimates, this imphes 



<c(||(/',r,r,5',5")iiL2 +e||(/,/',r,9,9')iiL^+^'lhlL2 ) 

21 

-I- fc 
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<C{e\\{fJ'J\g,g')\\H^ +e^u\\^, ). 
Substituting in (|7.7|) . yields 

<C7(||(/,/',5)IL2 +e||(/,/',r,g,5')llHi +£'hllL2j- 

Hence for e small, 

(7.10) e||n||^._^ <C(||(/,/^5)||z.^2^^+e||(/,/^/'^5,5')llHl,)• 
Plugging this estimate in (j7.8p 

(7.11) \\U'\\^, +\Hl^ +44l^ <C\\{f,fJ",g,g')\\H^+). 

£,S £,S e,S ^t" 

Hence, with (|7.9p . one has 

\\u"h2 < 

(7.12) 

C{ {f'J"J"',g',g")\\L^+e {f,f,f",g,g')\\Hi). 

Therefore, 

(7-13) \\U'\\^. + II^IL. < C\\f, f , f" ,g,g'\\^. 

e.d £,o £,o £,o 

The left hand side dominates 

and the right hand side is smaher than or equal to ||/||„2 + H^H ui • The estimate (17. 6p 

£,6 £,6 

follows. □ 
Knowing a bound for ||u||r2 , Proposition 16.31 immediately implies 

Proposition 7.3. There are constants C , > and 6q > and for s > 3 there is a 
constant Cg such that for e g]0, eo]; ^ G [0, (5o], / G ^e^^ > 9 ^ ^eS ^'^'^ ^ ^ s satisfying 
()5.8p and ()5.9p . one has 



(7.14) P\\H^_^^C^~'{\\f\\Hiy + \\9\\Hs_)+Cs\\U\\^s-i. 
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7.3 Proof of Proposition 17.11 

By Assumption I2.4r i). we may assume that there are hnear coordinates u = {ui,U2) € 
M"i X M"2 and h= lhi,h2) e x M"^, with n2 = rank 6*(n) such that 

(7.15) Kiu) = (, ° . , ° . 

\b2i{u) b22{u) 

and 622 ('u) is uniformly invertible on U^,. Introducing the new variable 

(7.16) U2 = U2 + Vui, V = {b^^)-%i{uNs), 
the equation b^u' — df^u = h has the form: 

a^^ui + a^'^U2 = h-\ , 

(7.17) 

V' V 1:22-/ -21 -22- J 

U2 — a ui — a U2 = ri2 



where 



-^■■=df*{% id)+^*(^' 0) 



Assumption I2.4l fii) implies that the left upper corner block a^^ is uniformly invertible. 
Solving the first equation for ui, we obtain the reduced nondegenerate ordinary differential 
equation 

i:22~/ , ^21^-11n-1-12 ~ -22- _ i , ^21,'-11n-1! 



6f n'2 + a^'ia'T a''u2 - a''u2 = h2 + a^'ia'T hi 



or 



(7.18) bu'2-du2 = h = 0{\hi\ + \h2\). 

Note that det J/* = det a^^ det d by standard block determinant identities, so that 
detd ~ det df^ by Assumption I2.4l fii). Moreover, as established in |MaZ4j . by Assump- 
tion [2?6] and the construction of the profile u^s we find that m := {b)~^a has the following 
properties: 

i) with m± denoting the end points values of m, there is > such that for all k : 

(7.19) \d!;{m{x) - m±)| < e/^+ig-^^l^l; 

ii) m{x) has a single simple eigenvalue of order e, dented by efi{x), and there is c > 
such that for all x and e the other eigenvalues A satisfy |Re A| > c; 

iii) the end point values of fi satisfy 

(7.20) fi->a /i+ < -a 
for some a > independent of e. 
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In the strictly parabolic case det 6^, 7^ 0, this follows by a lemma of Majda and Pego 
[MP]. 

At this point, we have reduced to the case 

(7.21) u'2-m{x)u2 = 0{\hi\ + \h2\), 

with m having the properties listed above. The important feature is that m' = O(e^) << e, 
the spectral gap between stable, unstable, and e-order subspaces of m. The conditions 
above imply that there is a matrix lo such that 

p := uj^^moj = blockdiaglp"*", e//,p— }, 

where the spectrum of p± lies in itRe A > c. Moreover, uj and p satisfies estimates similar 
to (I7.19p . The change of variables U2 = ojz reduces (I7.2ip to 

(7.22) z' -pz = uj'^Jz + 0(|/ii| + |/i2|). 

The equations (z^)' — p'^ z'^ = and {z^)' — p^ z^ = either by standard linear 
theory |Hej or by symmetrizer estimates as in |GMWZ] , admit unique solutions in weighted 
spaces, satisfying 

||e^l^lz±||i2 < C||e''l^'l/i±||i2, 

provided that 5 remains small, typically 5 < |Re p^\. 

The equation ZQ — £fiZQ = ho may be converted by the change of coordinates x — > x := ex 

to 

(7.23) dxZo - fj.{x)zQ = ho{x) = £'^ho{x/e), 
where zo{x) = zo{x/e) and /i(x) := ^{x/e). By (I7.19P 

|A(5) - Mil < Ce-^'*l 

with /i-t satisfying (j7.20p . This equation is underdetermined with index one, reflecting the 
translation-invariance of the underlying equations. However, the operator dx — /i has a 
bounded right inverse {dx — A)"^, as may be seen by adjoining an additional artificial 
constraint 

(7.24) 5o(0) = 

fixing the phase. This can be seen by solving explicitly the equation or applying the gap 
lemma of |MeZ3| to reduce the problem to two constant-coefficient equations on x ^ 0, with 
boundary conditions at z = 0. We obtain as a result that 

if 5 < minja,^}, which yields by rescaling the estimate 

lle^'^l^l zo\\l2 < Ce~^\\e'^\''\ho\\L2 
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Together with the (better) previous estimates, this gives existence and uniqueness for 
the equation 

z' — pz = h, Zq{Q) = 

with the estimate ||e^'^l^l2;||/,2 < Ce~^||e'^''l^l/i||^2. Because uj^'^oj' = O(e^), this imphes that 
for £ smah enough, the equation (|7.22|) with 2;o(0) = has a unique solution. Tracing back 
to the original variables n, the condition 20(0) = translates into a condition of the form 
4 • n(0) = 0. Therefore, the equation h^,u' — df^u = h has a unique solution such u that 
is ■ u{0) = 0, which satisfies 

||e"^l^lu||i2 < Ce-i||e^^'^'/i||i2 
for 6 and e small enough, finishing the proof of Proposition 17. 1[ 

Remark 7.4. The estimate of Proposition 17.11 may be recognized as somewhat similar to 
the estimates of Goodman [Goj in the time-evolutionary case. More precisely, the argument 
is a simplified version of the one used by Plaza and Zumbrun [PZ] to show time-evolutionary 
stability of general small-amplitude waves. 

Remark 7.5. The argument of Proposition 17.11 indicates that the estimate may be im- 
proved by factor e in transverse modes z±. However, we see no way to use this to improve 
the overall estimates on our iteration scheme. 



8 Existence for the linearized problem 

The desired estimates (I5.13P and (I5.14p are given by Propositions 17.21 and 17.31 It remains 
to prove existence for the linearized problem with phase condition 'u(O) • r(e) = 0. This we 
carry out using a vanishing viscosity argument. 

Fixing e, consider in place of C^U = F the family of modified equations 

(8.1) £^.''t/:=£^[/-r?(^^)=F:=Q, 4 • n(0) = 0. 
Differentiating the first equation yields 

(8.2) AU' - dQ{x)U - U" = (^'^ , 4 • u{{)) = 0. 
where dQ{x) denotes here the matrix dQ{uN s ■,v^{un s)) ■ 

8.1 Uniform estimates 

We first prove uniform a-priori estimates. We denote by the Schwartz space and for 
> 0, by the space of functions u such that e'^^^^^u G ,5^, with (x) = VT+~x^ as in 
([5TT]) . 
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Proposition 8.1. There are constants Eq > 0, 6q > and rjQ > 0, and for all s > 2 
a constant Cg, such that for e G]0,eo]! ^ £ [0,5o]> f] S]0,r7o], and U and F in ^^^(R), 
satisfying (|8.1|) 

(8.3) \\U\\^s <Cse~W\f\\HS+i + \\g\\j,s )• 

£.S e^o £,S 

Proof. The argument of Proposition 16.11 goes through essentiahy unchanged, with new rj 
terms providing additional favorable higher-derivative terms sufficient to absorb new higher- 
derivative errors coming from the Kawashima part. More precisely, consider again the 
change of variables v ^ v = v + pu, p = dvq~^duq{uj^s, v^,{u]\Js))■ Denoting U = {u, v) and 
U = P{u]\fs)U, (|8.ip is transformed to 

(8.4) AU' -QU - r]U" = 

with A, Q as in (16.60 . h given by (I6.10p and g now defined by 

9 = 9 + R2if + e'Cam + 7?(2p'n' + p"n). 

Thus we are led to equations of the form (16. lip with the additional term —r]U" in 
the left hand side. Using the symmetrizer S ()6.15p . one gains ^||C^"||^2 + -^11 1^' 1 1^2 in the 
minorization of Re {SF, U) and loses commutator terms which are dominated by 

v\\S"\\l^{\\U'\\1, + \\U\\l2\\U'\\l2)+ii\\K\\l^{\\U'\\l2 + \\U\\l2)\\U"\\l2, 

which can be absorbed by the left hand side yielding uniform estimates 

(8.5) V^\\U"\\l^ + WU'Wl^ + \\v\\l^ < C{\\f\\m + \\h\\m + \\~g\\m+e\\u\\L2). 
Going back to ()8.2p , this implies uniform estimates of the form 

(8.6) ^\\u"\^, +\\U'\\^, +\\v\\^, <C{\\{f,f\f",g,g')\\L. +e\\u\\^,). 

for 5 = 0, and next for 5 G [0, 5q\ with (5o > small, as in the proof of Proposition 16.11 

When commuting derivatives to the equation, the additional term rjd"^ brings no new 
term and the proof of Proposition 16.31 can be repeated without changes, yielding estimates 
of the form 

^\\dIu"\\l. +\\dlu%2 +\\d^ML\ 

£.6 £,6 £,6 

(8.7) <C\\d'AfJ',f",9,9')\\Ll, 

+ e''Ck{\\U'\\^k-i + e\\v\\^k-i + e\\u\\L2 ) . 

£,3 £.6 

Next, applying the Chapman-Enskog argument of Section [7] to the viscous system, we 
obtain in place of (j7.3p the equation 

(8.8) Ku' - df,u = f + 0{\v'\ + \g\ + l/'l) + e^O{\u\) + r^O{\u'\ + \U"\), 



f' + eh 
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where the final rj term commg from artificial viscosity is treated as a somxe. One applies 
Proposition 17. II to estimate e||ii||^2 ^ by the ^-norm of the right hand side, and continuing 

as in the proof of Proposition 17. 2^ the estimate (j7.13p is now replaced by 



^^■^^ < C{\\f,f',f",9,9%, + \\U"h2)). 



Therefore, for rj small, the new 0{rj) terms can be absorbed, and (j8.3p for s = 2 follows as 
before. The higher order estimates follow from (18. 7p . □ 

8.2 Existence 

We now prove existence and uniqueness for (|8.ip . First, recast the the problem as a first- 
order system 

(8.10) W -MA = r 

with 

U = \ V \ , T 



and 




/ ^11 Axi 

(8.11) A := 7?-M 

\vr^A2xAx\ - Q21 r]'^A2iAi2 - Q22 

Next, consider this as a transmission problem or a doubled boundary value problem on 
X ^ 0, with boundary condtitions given by the n + 2r matching conditions U{0~) = U{0'^) 
at x = together with the phase condition • u(0) = 0, that is n + 2r + 1 conditions in all: 

(8.12) U{0-)=U{0+), £e-u{0) = 0. 

Note that the coefficient matrix A converges exponentially to its endstates at ±00. 

Lemma 8.2. There is 9i > such that for e small enough , the matrices A± have no 
eigenvalue in the strip |Re z\ < e6o. 

Proof. The proof is parallel to the proof of the estimates. Dropping the ±, suppose that ir 
is an eigenvalue of A, or equivalently that there is a constant vector U ^ such that e^'^^U 
is a solution of of equations ()8.ip Thus 

Aiiu + A12V = irrju, 

(8.13) 

{irA -Q + T^r])U = 0. 
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Introduce once again the variable v = v + Q22 Q21U, so that the equation is transformed 

to 

Al^u + A12V = irrju, 

(8.14) . 

(irA-Q^ + T^r])U = 0. 

where A and Q now denote the end point values of the matrices defined at (j6.6p . Denoting 
by S and K the end point values of the symmetrizer and Kawashima's multipliers associated 
to A and Q, consider the multiplier 

S = \t\^S -irK - XS. 

Multiplying the second equation in (j8.14p by S and taking the real part of the scalar product 
with U yields 

IrpRe {KA - SQU, U) + X(SQU, U) + r]\r\'^{SU, U) 
<C(|ImT|(|r|2 + A))|C/|2 + C|r||QC/||[/| 

+ r?(|r|2|ImT|2 + |r|3 + A|r|2))|C/|2. 

Therefore, choosing appropriately A, for r/ and |Imr| sufficiently small, one has 

(8.15) (r/|r|^ + \t\^)\U\^ + \v\^ < C\linT\\u\^ 

In particular, |r| must be small if Imr is small. 

Prom the equation iTA2iu + A22V — Q22V + ^t^u = one deduces that 

V - iT{Q22)~^A2lU = 0{\t\ + r]\T\'^)\v\. 

Substituting in the first equation of (j8.14p , we obtain the Chapman-Enskog approximation 

{Al^ - iTK)u = 0{r]\T\ + |r| + r/|rp)|Imr|i))|n| 

where 6* denotes the end point value of the function p.9p . Therefore, 

(8.16) \{K)-'^A\^u - iTu\ < C|ImT|^|r||n| 

with arbitrarily small c > 0. We know from Assumption 12.61 that for e small, (h^)^^A\^ 
has a unique small eigenvalue, of order 0(e), real. Let us denote it by e^. Then we know 
that is bounded from below, see (j7.20p . Then ()8.16p implies that there is a constant C 
such that for |Imr| small enough, and thus |r| small, |ir — e//| < C|Imr|2 |r|. Therefore, 
/Imr + < ^e\fi\ if e is small enough. 

Summing up, we have proved that if e is small enough, A has at most one eigenvalue z 
in the strip |Re2: < e2|;u|, such that \z — < This implies the lemma. □ 
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Remark 8.3. The same reasoning can be applied to prove that A actually has a simple 
eigenvalue such that |z — < ^e|/i|. 

Proposition 8.4. There are constants Eq > 0, 5q > and rjQ > such that for e G]0,eo], 

5 G [0, 5o]) V £]0,r/o]) CLf^d F in ^^^(M), (|8.ip admits a unique solution U G ^(.^(M). 

Proof. Noting that the coefficient matrix A converges exponentially to A± at ±00, we may 
apply the conjugation lemma of [MeZlj to convert the equation (|8.1U|) by an asymptotically 
trivial change of coordinates U = T{x)Z to a constant-coefficient problems 

(8.17) Z'_- k_Z_ = Z'^- K+Z+ = F+, 

on {ibx > 0}, with n + 2r + 1 modified boundary conditions determined by the value of the 
transformation T at x = 0, where A± := A(ibcxD), and Z±[x) := Z{x) for itx > 0. 

By standard boundary- value theory (see, e.g., jHej ) . to prove existence and uniqueness 
in the Schwartz space for the problem (j8.10p on {x < 0} and {x > 0} with transmission 
conditions (|8.12p . it is sufficient to show that 

(i) the limiting coefficient matrices A± are hyperbolic, i.e., have no pure imaginary 
eigenvalues, 

(ii) the number of boundary conditions is equal to the number of stable (i.e., negative 
real part) eigenvalues of A_)_ plus the number of unstable eigenvalues (i.e., positive real part) 
of A_ , and 

(iii) there exists no nontrivial solution of the homogeneous equation f = 0, g = 0. 
Moreover, since the eigenvalues of A± are located in {|Rez| > 6ie, the conjugated form 

()8.17p of the equation show that if the source term / has an exponential decay e~^^^^^ 
at infinity, then the bounded solution also has the same exponential decay, provided that 

6 < 9i . Therefore, the three conditions above are also sufficient to prove existence and 
uniqueness in .5^£5 if e and 6 are small. 

Note that (i) is a consequence of Lemma 18. 2^ while (iii) follows from the estimate ()8.3p . 
To verify (ii), it is enough to establish the formulae 

dim5(A±) = r + dim5(Atf), 
^' ' dimZ^(A±) = r + dim^/(^*f), 

where A^f = dftf{u±) = An + Ai2dv*{u±) and S{M) and U{M) denote the stable and 
unstable subspaces of a matrix M. We note that ^4^^ = df^:{u±) are invertible, with 
dimensions of the stable subspace of A'^ and the unstable subspace of A^ summing to 
n -f 1, by Proposition 13.21 Thus, (j8.18p implies that 

dimcS(A+) + dimZY(A_) = 2r + d{mS{Al'[) + dim^^(A*j;) = 2r + n + l 

as claimed. 



28 



To establish (jS.lSp . introduce the variable v = v + Q22 and the variable corre- 

sponding to v' scaled by a factor rj^ ^ that is it; = 772^; + r]~2Q~^Q2i(Aiiu + A^v). After 
this change of variables, the matrix A it conjugated to A with 

(8.19) 7?5A= / . 

\0 -Q22 0/ VO(r?-^) 0{7]-^) A22) 

From (i), the matrix r/aA has no eigenvelue on the imaginary axis, and the number of 
eigenvalues in {Re A > 0} is independent of r/, and thus can be determined taking r/ to 
infinity. The limiting matrix has r eigenvalues in {Re A > 0}, r eigenvalues in {Re A < 0} 
and the eigenvalue with multiplicity n, since —Q22 has its spectrum in {Re A > 0}. The 
classical perturbation theory as in |MaZl] shows that for r/~2 small, rya A has n eigenvalues 
of order r]^^, close to the spectrum of A^^ with error 0(r/^^). Thus, for 7] > large, 772 A 
has r + dimS{A*^) eigenvalue in {Re A < 0}, proving (|8.18|) . 

The proof of the Proposition is now complete. □ 



8.3 Proof of Proposition 15.21 

Let (£*''')^ denote the inverse operator of Cl'^ defined by (18. ip . for sufficiently small rj > 0. 
The uniform bound (|8.3p . and weak compactness of the unit ball in H^, for F £ S^, we 
obtain existence of a weak solution U G of 

(8.20) ClU = F:= Q , 4 • ^^(0) = 0, 

along some weakly convergent subsequence. Proposition 17.21 implies uniqueness in for 
this problem, therefore the full family converges, giving sense to the definition 

(8.21) (rjt = lim(£:''')t 
acting from to H^. 

For F G Si;S, the uniform bounds (18. 3p imply that the limit (£^)^[7 G ^ and satisfies 
same estimate. By density, the operator (>C^)^ extends to / G and g G with 

(£J)tF G HI,. 

The sharp bound (|5.13p and (j5.14p now follow immediately from Propositions 17.21 and 
17.31 The proof of Proposition 15.21 is now complete. 

Remark 8.5. We have used freely the finite-dimensionality of v in our proof of linearized 
existence. However, as promised, it plays no role in the final linearized bounds. Thus, our 
result may be used together with discretization (Galerkin approximation) of v to obtain 
results also in the case that v is infinite-dimensional, as we do for the Boltzmann equations 
in [MeZ2j . 
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9 Application to spectral stability 



Proof of Corollary 14-^ In |MaZ3j . under the same structural conditions assumed here, it 
was shown that smaU-amphtude profiles of general quasilinear relaxation systems are spec- 
trally stable, provided that 



(9.1) 
and 

(9.2) 



\U' 



<c\u. 



El 



- + sgn(?7)i?o 



\U"{x)\ <C\U+-U-\ \U'{x)\, 



<C\U+- U^\, 



Rn 



r{uo) 
^dv^.{Uo)r{uo)^ 

where r(uo) as defined in Theorem 14.11 is the eigenvector of d/* base point Uq in the 
principal direction of the shock. From the bounds of Theorem 14. H we immediately verify 
these conditions, giving the result. □ 
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